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$1. Introduction. 


N two former papers t, which will here be referred to 
as I. and II., the general theory of the diffraction of 
electromagnetic waves from an electrical point source 
round a finitely conducting sphere was considered in detail. 
For the application to radio communication general 
formule for the field at any distance from the transmitter 
were derived. Special practical formule were given for 
two limiting cases: (a) total reflexion corresponding to 
an infinitely conducting sphere ( |k, | >œ) and (b) 
maximum absorption corresponding to 


sie). ae 
zA klk 2—1 
where the same notation as in the former papers is used. 


It is the purpose of the present paper, while making 
full use of our former results, to derive approximations 


<<x<I, 


* Communicated by the Authors. 
t Balth. van der Pol and H. Bremmer, Phil. Mag. xxiv. pp. 141 and 
825 (1937). 
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suitable for numerical application for any arbitrary con- 
ductivity c and dielectric constant «e of the earth and for 
any frequency. In order to simplify the problem we 
assume both sender and receiver to be situated on the 
earth’s surface (h;=h,=—0). While adhering to our 
former notations, the present investigation, moreover, 
will enable us to compare our numerical results with those 
of other writers. 

The basis of our present investigation will, therefore, 
be the series I. (15) for the Hertzian vector M,e, but for 
b=a, which thus reduces to 


24/27 Etim sa nil2 
VJ be „éS, (1) 


Hyo ~l ~ ——e o anes 
or 8" (ka) sing "=o (Sn 


where (see I. (6)) 


1 
a. 


Nee) = [5 e ogeta] 


d 
~ [p E Be} Teme’ ® 


In these formule n, represents the zeros with positive 
imaginary part of N,, ..;,2(z, 2’) considered as a function of n. 
Our problem, therefore, consists of a further numerical 
determination of the position of these zeros together 
with their corresponding residues. 


§2. The Zeros of Nn_1,(z, 2’) considered 
as @ Function of n. 


In order to obtain a clear picture of the position of the 
individual zeros of N,_,,., we first decompose the out- 
going wave %4? a(z) occurring in Nai, as the complex 
sum of two standing waves p-n- and p12 according to 
the known formula 


ge ceed ai 
3 ~ 3 
aa) o a Pte nl 


Ld e e (3) 
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Substitution of (3) in N,_,,.=0 leads to the following 
linear relation in e”™ : 


2b 1,2(2’) $ 7 feab, 1;2(2)} —2 bn) (2) = {z'n 12(2')} 
z'aniplz’) ë Seni} Pn ark) rA (paale) 


(e=ka; z? =k), . (4) 


envi 


which constitutes an equation in n. Now, in order to 
show clearly the position of the different roots n=n,, we 
write (4) in the equivalent logarithmic form 


= — Pn,— 122) ” , \ x = 
(n, —l}mi=log { haapt e R” m-i p H2mis, (5) 


where the logarithmic term is meant to represent the 
principal value, whereas the multiplicity of the logarithm 
is accounted for by the term 2ris (s=0, 1, 2...), and 
where in (5) we introduced a new reflexion coefficient 


= E z log TAO 


+ E 7 log ele) i, 


z=kya 


=kea 


R"(n) = ——;--——— - 


(6) 
z=ha 


_ [ 3 log TADIN 


It is worth noticing that this reflexion coefficient built 
up of % functions only, and therefore suggesting standing 
waves, is different from the R,, of II. (68 a), in which 
¢ functions (i.e., progressing waves) occurred only. 
It also differs from R, of J. (2) which was of a mixed type. 
Tn a later paper we hope to consider in more detail the 
physical character of these different reflexion coefficients 
which are closely related to the dual d’Alembert and 
Bernoulli solutions of the wave equation in one or more 
dimensions. 

The expression (5) is still rigorous, but two approxima- 
tions for the Bessel functions, already considered before, 
are available : I. (16) and I. (17) where, as shown in I., the 


3H2 
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expression (17) leads to more accurate results. Considering 
the less accurate case I. (16) first and making use of 


Pn —12(%) = HE l) Saale) 
pnl) = Herml (2) pe" L? peh 
one finds for I,n>0, |n | and |z|>>1, 


y 


d n* 
de log {zni} ~ v =l, 


2 
Z iog {2b n—12(2)}~— vi oad, 


22 


bin—1/2(2) Pe ais i" Vina 
eT a A 4 $ 
Y_n—12(2) 


Substitution of these expressions in (5), and calling 
n,=G@A,, leads to 
kis ,---- du i 
ian. | vi—u z + g -log R(As)=(s+4)m, - (7) 
1 
where 
RA \= keV Age ktkt Ag? Ie? E ig(3' —9") 

PO beg? gk tke Vie — hy? tg +Y 
and the spherical reflexion coefficient R” thus degenerates 
to the Fresnel plane reflexion coefficient as occurring in 
the plane Sommerfeld problem. 5’ and 3” represent the 
angles of incidence and refraction according to the 
relations A=k, sin 3’=k, sin 3”. 

Our result (7) has the character of a phase-integral 
equation and is the same as that derived by T. L. 
Eckersley *. However, it should be noted that the 
term 7/4 in the right-hand member of (7) appears here in 
a stringent way, whereas in Kckersley’s considerations, 
which leave an absolute phase constant undetermined, 
this term is borrowed from Watson’s analysis. 

Returning to (7), it may be remarked that for practical 
purposes the roots %, ny, ng of low order s are of domin- 
ating importance, because they are the roots with the 
smallest imaginary part, representing in the final ex- 
pression the terms with the smallest attenuation. Thus 


* T. L. Eckersley, Proc. Roy. Soc. London A, exxxvi. p. 499 (1932), 
see pp. 520-523, where Eckersley’s variable s,should be replaced by 
s—1 for comparison with our notations. 
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considering (7) again as an equation in n, we obtain the 
approximation 

Ns =2-+-2""7,, 
where 7, satisfies 


—4(—27,)?? +arc toV —27,)==(8s-+3)0,. . (8) 
where 
; dR ika?k; 
ENN EOE sie is a 9 
22"? ds! Pe 23V keke —1 i 
-and 


0<Refaretgl <r, 
—n/2<arg V—2r<0. 
Our present expression (8) fully confirms the two 
limiting cases of total reflexion and maximum absorption 


considered in I. and II., for, if in (8) we take |k, | >a 
and therefore | 5 | +00 , it degenerates to 


t~ h{8n(s +h) Prem, 


which is our old I. (40a), whereas for |è | <<1, (8) 
degenerates to 
te~ H 3n(8+2)}2 Feit 5, 


which is the other case already considered in IT. (page 853). 


The phase integral equation (7) and its approximation 
(8) for the zeros of low order s were derived with the aid 
of the approximations I. (16) for the Bessel functions. 
As already remarked, we gave in I. also another approxima- 
tion for Bessel functions, viz., I. (17), which is considerably 
more accurate, especially in the neighbourhood of n=z, 
and therefore gives us a better approximation for the low 
order zeros n, of Naija The numerical difference of these 
two approximations was already fully considered in I., 
page 163. The application of this better approximation 
to the Bessel functions occurring in the function N,_1/, 
leads to 
Bin V = 27 HY3{4( —27)??} l 


V ZET Hag <, (10 
ghi8 HYj3{4(~27) 37} 24138 ( b) 


Nica 


whereas the approximation I. (16) yields 


a V — 2r mn. (—2r)32 1 
Naza ~- ag 9 1t g } ps: (100) 


Z 
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Hence we obtain as equations determining the zeros n 
e- mH G (—2r, je 2} 1 
HGR RP TT SV 
which equation we will call the Hankel approximation 
and which may be compared with the following inversion 


of (8), 
m (— 2r)? n l N 


which we will call the tang. approximation. In order to 
obtain a clear distinction between expressions based on 
the tang. and Hankel approximation we shall also later on 
number the formule based on the first one with an a 
and those based on the second one with a b. It may be 
remarked that the older work of T. L. Eckersley * and 
also the more recent investigations of B. Wwedensky + 
are based on the fang. approximation. This was also 
the case in I. and II., but only when | |<<1, where, 
as will be explained further on (see (14)), the error in 7 is 
small and less than 1 per cent. 

Having thus treated the general aspect of the approxi- 
mations, we will in the next paragraph elaborate the 
numerical behaviour of the roots of (11). 


(11) 


$3. The. Numerical Values of the First 
Few Zeros of Ny-yj9(2, 2’). 


The behaviour of (11) is completely determined by the 
numerical value of the 8 of (9). This variable has in its 
denominator the great factor 2/8=(k, a)", and |8| may 
therefore in practice obtain values as low as 0-005. On 
the other hand, k/k; may reach such big values that it 
more than compensates the factor z—1/*, so that then | 6 | 
may be as high as 1000. We shall, therefore, have to 
evaluate the zeros of (11) for a great range of values of 6. 
As, moreover, we shall derive the zeros ns of N,_,,. for 
finite values of 5 from the two limiting cases $= œ and 
5=0, a very exact knowledge of the values of n,=2-+2"87, 
for these limiting cases is necessary. 


* T. L. Eckersley, loc. cit. 
+ B. Wwedensky, Techn. Phys. U.S.S.R. ii. p. 624 (1935) ; iii. p. 915 
1936); and iv. p. 579 (1937). 


Waves over a Finitely Conducting Spherical Earth. 823 


It follows from (11) that the first case, corresponding 
to total reflexion, leads to 


(27)? ==(3r)8'2(s+f)e-*,. . . (12a) 
or H{}{4(—27,)?}=0, which is equivalent to 
Jall — 27,)8?e™} —I_oia{3(—27,)87e™} =0, (126) 


whose roots we shall call 7, .., whereas the second case, 
corresponding to maximum absorption, leads to 


(—27,)"?=(37)'%(stRe-™, . . . (18a) 


or HR 


(){4(—27,)3?}-=-0, which is equivalent to 


Ti sf 2r in) + T_yofG(—27)"e}=0, (130) 


with the roots 7, o, say. 
A careful numerical evaluation gave us :— 


Tang. approximation. Hankel approximation. 
To, a = 0°885.e°"!8 To, » = 0:808. et 
Ti =2589. 673 Ti, o =2 577 ei 
See | To; o =3'830 e™3 To, m = 3'824 eT 3 
T3, o =4°895.68™3 T3, o =4'892 eim? 
T4, s= 0853. T3 T4, o = 5'851. e73 
T5, , =6°739.¢im/3 T5, o = 6737 ei 
a 
í ro p= 1842. e878 To 0= 1856.73 
Ti, o= 3° 240.68 71, 0==3 245. e™3 
jg o= 4379. eT To, o= 4'382 e878 
T3, o= 5384. eiT T3, o= 5'386. ei™/3 
T4, o=6'304. e773 74, o=6'305 , ei"!3 
T5, o= T 161.3 T5, o= 1161.68 


Now, in order to deduce the 7,’s for any finite 5 from 
these limiting values, we make use of the Riccati equations 
for 7 as function of 8, which can be derived from (10), 
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using the known properties of the Hankel and the tang. 
functions respectively, viz. : 


a 23% Pe ee . . . (14a) 
dr 2r 
mS EE E O . . (146) 
dr 


(14 b) can also ea be derived for |z | >>1 from 


Nn-ipe™ = É 1 log pape — 
with the a of 


gh a5 = 0, 


TERO she RED Lat @)=0 


It also follows from a comparison of (14a) and (14 b) 
that the error by using the tang. approximation is small 
only when |è | <<2 |r, |. 

Now that the limiting values 7, ,, and 7, are known, 
we can with the aid of (14 b) derive two series of asymptotic 
character, expressing 7, as a function of 5-! and of ô. 
These are easily found to be 


1 1 1 3 1 
—- Sali ele 
275, að BT, .°d? 127, Be ( T ETs, 2) 83 


o l ( E O. 
B27, o 3r at | Er, 98) 81 

i 21 35 \ 1 

~ Ae ETa + 3275, & 


o 1L 29 Bl f 
4875, o? {i5 hae a a he ee a) get 


2 l 
Te = Ts, poo 37% 50° + 2 t t p Ts otin, 006 


1 58 
Tin (5487, 08" + Tg Ts, Obon 
5b) 


Fortunately the numerical values derived from the first 
series for | 527 |> 4 and from the second one for | 5’7| <4 
yield a good fit even for complex values of ò. As, further, 
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the imaginary parts of the 7,’s will determine the attenua- 
tion of the different terms in our final expression for Tow 
we add the corresponding series (176) (186) for the 
imaginary part of the more important rọ and 7,. In 
these series we have also introduced 8 by its modulus and 
argument as 


8K.) (16) 


where the argument is taken as ?7—y, because then the 
case often encountered in practice, where the conduction 
current is great in comparison with the dielectric current, 
we have y=0. 

The two series referred to above are :— 


0-6189 sin (4% —15°) , 0°2370 cos (2%) 
K K2 
0°0539 sin (34+15°) , 0:00273 sin (444,—60°) 
K + n am 
900211 sin {(54+45°) | 000142 sin (6% —30°) 
KS hi Ke as 
Imz)=1-6070--sin (Y +45°)K— 1:237 sin (3%-+ 75°) K3 
+4 sin (44)K4+2°755 sin (54—75°)K$ 
—2°887 sin (6% —150°)K®— 6:589 sin (76+- 135°) K* 
— 13°32 sin (8—120°)K®..., 


Imr,=0°6996-+ dy 


(176) 
and 
0:1940 sin (—15°) , 0:0073 cos (2%) 
Im7,=2-2317-+ k oe p 
0-0120 sin (3% =+15°) _ 0°00160 sin (4%— 60°) 
D eet 7 
0-00124 sin (54+45°) _ 0-000324 sin (6/—30°) 
pe ke ee 


Im7,==2°810—sin (++ 45°)K — 2163 sin (3y +75°)K? 
+4 sin (4~)K4+8-422 sin (54 —75°)K5 
—5:047 sin (6—150°)K®— 38°32 sin (7+ 185°)K? 
— 40°71 sin (8¢—120°)K8. .. 
(18b) 
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These series are graphically represented in fig. 1, where 
the abscisse represent K= |8 | and the ordinates 
Imr, and Imr, the parameter being the phase ¢ of è 
as defined in (16). As % varies between 0 and 45°, we 


Fig. 1. 
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Imro and Imz, as a function of | 6 | =K. 


have chosen for this parameter the values 0, 15°, 30°, 
and 45° only. 
From this figure it may be noted, as also follows from 
(17), that for 4=0, Imr, shows a minimum as function of 
|è |. This fact, as we shall see later on, has important 
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numerical consequences in the theory of radio propagation. 
However, it should be observed that this minimum dis- 
appears when > 15°. 


Fig. 2. 
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We are now prepared to compare our result with a 
similar one obtained for Imro only by Wwedensky *. 


* B. Wwedensky, loc. cit. ii. pp. 633-635 (1935). 
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and we therefore replot the lower part of fig. 1 for J=0, 
but introduce the new, practical variables 


B= V21dy,. IMT = 34°20. Imry, 
and 
e = ates., 10-5. ] 8 | = 4415. 10-7.] 8 |. 


Ckm, 


As, moreover, as stated above, Wwedensky used the 
tang. approximation, we also represent this approximation 
in fig. 2. It clearly follows from this figure that the two 
approximations practically coincide for o}? „A%6 <2, but 
that an appreciable difference remains when o!!?,_,Aj/§ >5. 
Again, the curves given by Wwedensky and Eckersley do 
not show the minimum obtained above, to which we will 
refer later on. 


$4. The Final Formule for the Hertzian Vector Iio. 


Returning to our basic formula (1) for I, the last 
necessary step is to investigate 
(2a) =A ($) 
On N= Ng glis Or THT, 
This is easily done by differentiating (10 a) and (10 b) and 
applying known differential relations for the tang. and 


the H functions respectively together with (11). One 
thus obtains 


Bea) 
n= Ne 


i 1 4 
PEE ee PE E A 
an gpl ~2te+ 57, t =) y 30) 


Nna lji 1 
Substituting these results in (1) and introducing the 
variable y=2! 39=%/ 2r. anc already considered in I. 
and II., we obtain as our final formule 


@itsx 


Ilio ~p. V 2riy. a 


f Te I T (20a) 
ats Sop = RG 
3 257, ò? 
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eiT 8x 


Toy ~ 21 Tr . NE. $ P E A 


s=0 (< l 
(el | 


laz 


(20b) 


(20a) is essentially the equivalent of Wwedensky’s 
final expression *. However, apart from the fact that 
the 7,’s of our expression (20 b) are defined by (11 b), and 
not by (11 @) as is the case for (20 æ), for the reasons given 
above we consider (206) as the proper solution, this being 
based on the Hankel approximation, which is accurate 
up to terms of the order z-?°. 

A verification of (20 b) can be obtained by considering 
the two limiting cases (a) total reflexion ->œ and (b} 
maximum absorption 5-0, fully confirming our former 
expressions IT. (100). 

It ean also be proved that, just as in I. ($5), the general 
rigorous formula for the plane boundary could be con- 
sidered as a limiting case for a—>oo of the general rigorous 
expression I (1a) of the spherical problem, so can the 
well-known approximate p formula of Sommerfeld + be 
derived as similar limiting case (for a->o ) of our final 
approximate formula (206). We hope to return to this 
question in a later note. 


$5. Numerical Results. 


Now that we have obtained our final expression (20 b} 
in a form readily adapted for numerical evaluation, it 
was a relatively simple matter to construct practical 
curves for the field at the earth’s surface of an emitter 
also placed on the ground. For an emitter radiating 
P kilowatt we have for the field in practical units 

__ 310° 


E Tos. 
Sum a cw. 
lad ee | ee kw. |2 TL. rl 


With these units we have calculated the figs. 3, 4, 5, and 6 
showing curves representing the field due to a sender 
radiating one kilowatt as a function of distance. Figs. 3 
and 4 refer to a conductivity o=—410-" emu. and 


* B. Wwedensky, loc. cit. ii. p. 631, formula (6, 6), where it should be 
noted that Wwedensky uses instead of ô : M = 1/(31K). 

+ See, e.g., Balth. van der Pol, Zeits. für Hochfreq. Techn. xxxvii. 
p. 152 (1931). 
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dielectric constant e=80, corresponding to sea-water, 
whereas figs. 5 and 6 were calculated for c=10-8, e=4 
corresponding to average soil. The different curves 
in each of these figures cover the frequency range from 
15 ke./s. to 300 mc./s, (A=20 kilometres to 1 metre). Anew 
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Field strength of the ground wave of a transmitter of 1 kw. for 
propagation over sea-water for distances to 2000 km. 


interesting fact may be read off from these curves viz., 
that for long distances long waves are slightly better 
propagated over average soil than over sea-water, e. g., 
at D=2000 km. the field of a wave A=20 km. is 4 per cent. 
stronger over average soil than over sea-water. This 
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is a consequence of the minimum of Imr, as function of 

|è | (fig. 2) already referred to above. It follows from 
fig. 2, that for a given wave-length, at least as long as the 
dielectric current is small compared with the conduction 
current, always a conductivity can be found for which 
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Field strength of the ground wave of a transmitter of 1 kw. for 
propagation over sea-water for distances to 200 km. 


at long distances the propagation conditions are optimum. 
As this minimum of Imro occurs near 15.10-', this opti- 
mum conductivity is given by 
eee 0:225 

e.m.u. Nye 


„1071, 
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Finally,it may be remarked that in constructing these 
curves the use of three terms only of the series (20 b) 
gives at shorter distances already a good transition to 
the values calculated with the aid of the p formula for the 


Fig. 5. 
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Field strength of the ground wave of a transmitter of | kw. for 
propagation over average soil for distances to 2000 km. 


plane problem. Hence, the curves could be extended 
to short distances from the sender. 
Summary. 


An extension is given to the applications in the radio 
case of a general theory treated in two former papers 
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on the problem of diffraction of electro-magnetic waves 
round an absorbing sphere. Asymptotic series have been 
derived from which the numerical constants occurring 
in the final formula may be calculated for the field strength 
of a sender emitting waves over a homogeneous sphere 


Fig. 6. 
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Field strength of the ground wave of a transmitter of 1 kw. 
for propagation over average soil for distances to 200 km. 


with arbitrary electrical conductivity and dielectric 
constant ; the sender and receiver are hereby supposed 
to be situated on the earth. The general results are 
compared with those of Eckersley and Wwedensky, 
treating the same problem. Finally, curves are given 
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showing the field strength as a function of the distance 
for propagation over sea-water and over average soil. 


POSTSCRIPT. 


Hereby we call attention to the following errata in our 
former papers :— 


(1) Page 151, in formula (18) read I instead of I,. 
(2) Page 164, line 6 from the top, read z753 instead 
of 2-18, 
(3) Page 171, in formula (55) read z’/z instead of z/z’. 
(4) Page 835, line 2 from the bottom, read I° , , instead 
of Ti. 
(5) Page 838, at the end of the first line of formula (81 a) 
read Tox44 instead of 7ox4. 
(6) Page 838, in formula (81 b) read 
2K+4 2K+4 
x instead of X. 
v=3 A=3 
(7) Page 854, in the denominator of the last line of 
formula (102), read 7’? instead of +’22, 
(8) Page 860, in the second line of formula (109), read 
P` (cos 3) instead of P} (cos 3). 
(9) Page 861, line 7 from the top, read and instead of en. 


(10) Page 863, in the formula of line 10 from the top, 
read 9/128 instead of g/128. 


Eindhoven. 
20th January, 1938. 


LXXVI. The Evaporation of Atoms, Ions and Electrons 
from Tungsten. By A. L. REIMANN, Ph.D., D.Sc.* 
(Communication from the Research Staff of the M.O. 
Valve Co., Ltd., at the G.E.C. Research Laboratories, 
Wembley, England.) 


Introduction. 


HE Saha-Langmuir-Kingdon theory ‘© of the equi- 
librium between neutral atoms, positive ions and 
electrons in a gaseous mixture at absolute temperature 


* Communicated by C. C. Paterson. 


